Time-dependent density functional theory on a lattice 



(N 

o 

(N 
C 



M. Farzanehpour^Q and I. V. Tokatly 1 ' 2 -Q 

1 Nano-Bio Spectroscopy group and ETSF Scientific Development Centre, 
Departamento de Fisica de Materiales, Universidad del Pais Vasco UPV/EHU, E-20018 San Sebastian, Spain 
2 IKERBASQUE, Basque Foundation for Science, E-48011 Bilbao, Spain 
(Dated: March 8, 2013) 

A time-dependent density functional theory (TDDFT) for a quantum many-body system on a 
lattice is formulated rigorously. We prove the uniqueness of the density-to-potential mapping and 
demonstrate that a given density is u-representable if the initial many-body state and the density 
satisfy certain well defined conditions. In particular, we show that for a system evolving from its 
ground state any density with a continuous second time derivative is w-representable and therefore 
the lattice TDDFT is guaranteed to exist. The TDDFT existence and uniqueness theorem is valid 
for any connected lattice, independently of its size, geometry and/or spatial dimensionality. The 
general statements of the existence theorem are illustrated on a pedagogical exactly solvable example 
which displays all details and subtleties of the proof in a transparent form. In conclusion we briefly 
discuss remaining open problems and directions for a future research. 
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I. INTRODUCTION 

Time-dependent density functional theory (TDDFT) 
is rapidly becoming one of the most popular methods for 
modeling dynamics of realistic quantum many-particle 
systems. Nowadays TDDFT is routinely applied to dy- 
namical problems in condensed matter, atomic, and nu- 
clear physics, in quantum chemistry, and in nanoscience. 
Only within last few years several journals published spe- 
cial issues on TDDFT— ~— [see also a recent review, Ref. H[ 
and two TDDFT textbooks came out earlier this year 5,6 . 

The reasons for the popularity of TDDFT are obvi- 
ously the same as those for the ground state DFT. It 
reduces calculations of the density in a complicated in- 
teracting many-body system to solving a set of Hartrcc- 
type equations for a reference system of noninteracting 
Kohn-Sham (KS) particles, which simplifies the problem 
dramatically. A possibility of such a reduction rests on 
two fundamental mathematical statements: (i) a one-to- 
one mapping between the density and the external po- 
tential, meaning that two different (modulo a constant) 
potentials cannot produce the same density; and (ii) a 
-y-representability of the density, i. e. the existence of a 
potential that generates a given density. The first state- 
ment, known as the mapping theorem, guarantees that 
a many-body wave function and thus any observable are 
unique functionals of the density. The w-representability 
is required to justify the KS construction in (TD)DFT. 

In the static DFT a proof of the mapping theorem 
was presented in a seminal paper by Hohenberg and 
Kohn 7 , while a general solution of a more tricky v- 
representability problem appeared much later and only 
for lattice systems^ which is probably sufficient in most 
practical cases. Proving corresponding theorems for 
TDDFT turned out to be much more difficult, mostly 
because of the absence of a minimum principle in dy- 
namics. Only in 1985 Runge and Gross (RG) succeeded 
to find a sufficiently general proof of the TDDFT map- 



ping theorem for a class of analytic in time (i- analytic) 
potentials^ .. An attempt to attack the time-dependent 
w-representability problem has been performed in Ref. [ll| 
by assuming i-analyticity both for potentials and for al- 
lowed densities. Under this restriction a formal power se- 
ries for the potential can be uniquely reconstructed from 
a given Taylor expansion of the density. Unfortunately 
the convergence of that series is unproved up to now, and 
thus a complete solution of the i>-representability prob- 
lem within the series expansion technique is still lacking. 
We note that the issues of i-analyticity and a uniform 
convergence of power series in quantum dynamics are not 
as exotic as it may appear on the first sighti2r— . Despite 
a number of indication s 15 : 16 that i-analyticity was not a 
fundamental limitation of the theory, a question of a more 
general and clean justification of TDDFT remained open 
for many years. 

Very recently it has been recognize d 13 : 1 7 ~ — that the ex- 
istence of all TDDFT-type theories is equivalent to the 
solvability of a certain universal nonlinear many-body 
problem which determines the potential and the many- 
body wave function in terms of a given basic observable. 
Mathematically this universal problem can be posed in 
two different forms. The first possibility is to view it 
as a Cauchy problem for a special nonlinear Schrddinger 
equation (NLSE) 13 : 17 ~ — . The uniqueness and the exis- 
tence of solutions to this NLSE are equivalent, respec- 
tively, to the mapping and the u-representability theo- 
rems in TDDFT. Alternatively it can be formulated as 
a fixed point problem for a certain nonlinear map in the 
space of potentials 2 ^. First applications of the above two 
formulations appeared almost simultaneously in the last 
yea r 20 ' 21 . 

The fixed point approach has been used in Ref. [2l| 
to prove the existence of the "classical" TDDFT in its 
original RG form. In this work the i-analyticity require- 
ment was completely relaxed and effectively replaced by 
a more physical and plausible assumption - a boundness 
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of a certain generalized response function related to a 
stress-density correlator. 

A rigorous formulation of a time-dependent current 
density functional theory (TDCDFT) on a lattice was 
presented in Ref. within NLSE formulation of the 
problem. The lattice TDCDFT turned out to the 
first and, in fact, the simplest example of a TDDFT- 
type theory for which both the mapping and the v- 
representability theorems have been proved without any 
unjustified assumption. 

In the present paper we further develop and extend 
the lattice NLSE technique of Ref. |2CJ to address a long 
standing problem of TDDFT for lattice many-body sys- 
tems. It may appear surprising, but even the RG map- 
ping theorem for i-analytic potentials is currently absent 
for the lattice TDDFT. This should be contrasted to the 
lattice TDCDFT where the standard power series argu- 
mentation can be easily adapted 2 ^. Many discussions 
of mathematical and conceptual difficulties of the lattice 
TDDFT can be found in the literature 2 ^— . In this work 
we show how to overcome these problems using NLSE 
formalism. We prove the uniqueness and existence theo- 
rems for the lattice TDDFT and analyze the conditions 
one has to impose on the initial state and the density 
to guarantee v-representability. In particular we demon- 
strate that practically any properly normalized density 
is locally v-representable provided the dynamics starts 
from the ground state. These results put applications of 
TDDFT to various lattice models^— on a firm ground, 
and shed new light on the general mathematical structure 
of TDDFT. 

The structure of the paper is the following. In Sec. II 
we present a general formulation of the lattice many-body 
theory and derive a lattice analog of the force balance 
equation that plays a key role in the NLSE formalism 
for TDDFT. Section II is the central part of the present 
paper. Section III A starts with a brief review of our ap- 
proach to TDDFT in a more familiar continuum case. 
Then we derive the corresponding NLSE for a lattice 
theory, and, finally, formulate and prove the general exis- 
tence a uniqueness theorem on the lattice TDDFT. Sev- 
eral important aspects of the basic theorem for a generic 
initial state are discussed in Sec. IIIB. Section IV presents 
an explicit illustration of the general existence theorem 
for a simple exactly solvable model - one particle on a 
two-site lattice. In Sec. V we consider a practically im- 
portant case of a system evolving from its ground state. 
The main outcome of this section is a u-representability 
theorem for the initial ground state. In Conclusion we 
review our results and discuss open questions and di- 
rections for further research. In Appendix we derive an 
explicit NLSE for another illustrative model describing 
dynamics of two interacting particles on two sites (a Hub- 
bard dimer). 



II. PRELIMINARIES: MANY-BODY PROBLEM 
ON A LATTICE 

In this work we are considering quantum dynamics of 
N interacting particles on a lattice which consists of a 
finite, but arbitrary large number M of sites. The state 
of the system at time t is characterized by a many-body 
wave function ip(r, r 2 , rjv; t), where coordinates r$ of 
particles (i = 1,2,..,N) take values on the lattice sites. 
The dynamics driven by an external scalar potential 
v(r;t) is described by the following discrete Schrodinger 
equation 

N 

id t ip(ri, ...,r N ;t) = - ^ ^T rijXj •••!*) 

i—l x, 
JV 

+ ^2'"(T j ;t)i)(Ti,...,r N ;t) 

3=1 

+ ^2w Ti ,r j il>(r 1 ,...,r N ;t), (1) 

j>i 

where real coefficients T T y = T r > r correspond to the 
rate of hopping from site r to the site r' (we assume that 
T r r = 0), and w r ^ T > is a potential of a pairwise particle- 
particle interaction. To cover various possible physical 
applications we do not specify the geometry of the lattice 
and do not assume as usual that the interaction depends 
only on the distance between particles. For example, 
the latter is important in a typical transport setup with 
noninteracting or weakly interacting leads connected to 
a strongly interacting central regio n 29 ' 31 " — . 

Equation ([T]) (formally it corresponds to a system of 
M linear ordinary differential equations) determines 
the wave function as a unique functional of the external 
potential and a given initial state, 

ip(r 1 ,r 2 ,...,r N ;t ) = ^ Q (r 1 ,r 2 , ...,r N ). (2) 

The key object of DFT is the density of particles 
n(r;t), which in the present context means the number 
of particles on a given site 

n(r;t)=N £ |^(r, r 2 , r N ; t)\\ (3) 

r 2 ,...,riv 

where we assumed that the particles are identical. By 
taking the time derivative of the definition and us- 
ing the Schrodinger equation ([T]), we find the following 
equation of motion for the density 

n(r; t) = i^[T r , r ^(r, r'; t) - T r ,, r p(r', r; t)}, (4) 

where n = dtn, and p(r,r';i) is a density matrix (or a 
"link density") on the [r, r']-link, 

p(r,r';t)=N ^ tp*(r, r 2 , r^; f)V>(r', r 2 , r N ; t). 

r 2 ,—,TM 

(5) 
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Obviously Eq. ((4]) is a lattice version of the continuity 
equation. Since in the left hand side of Eq. (U) we have 
the time derivative of the on-site number of particles, 
the right hand side should be identified with a sum of 
currents flowing along links attached to the site. Indeed 
introducing a link current from site r to site r' as follows 

J(r,r')=2Im[T r , I , /9 (r,r';t)], (6) 

we can rewrite the continuity equation (|4]) in more famil- 
iar way 

n(r;i)=-E J ( r < r ')- (7) 

r' 

This equation shows that the decrease rate of the density 
on each site equals to the sum of all outgoing currents. 
Equation ([7]) can be also viewed as an integral of the 
usual differential continuity equation over a small volume 
element surrounding the site r. 

Now we introduce another equation of primary impor- 
tance for the lattice TDDFT. This is a lattice analog of 
a divergence of the local force balance equation. It can 
be derived by differentiating the continuity equation Q 
with respect to time, and using Eq. (fTJ) to transform the 
derivative of the right hand side. After straightforward 
calculations the lattice force balance equation reduces to 
the following form 

n(r; t) = 2ReJ2 T r,*p(r, r'; t)[v(r'; t) - v(r; t)] - 9 (r; t). 

r' 

(8) 

Here q(r; t) stands for a " lattice divergence" of an internal 
stress force, 

g(r; t) = 2Re ^ T r , r -{p 2 (r, r", r'; t)(w r ,v» - «v, r ») 

r ; ,r /; 

+ lT r ,y,p(r,r";t)-T T y,p(r',r";t)}}, (9) 

where ^(r, r", r'; t) in the right hand side is the two body 
density matrix, 

p 2 (v,v",r';t)=N(N-l) £ </>*(r, r", r N ; t) 

r 3 ,...,r N 

x V(r',r",...,r^;t). (10) 

A special role of Eq. © for TDDFT follows from the 
fact that it explicitly relates the potential v(r;t) to the 
density n(r; t) and the instantaneous many-body state 

The force balance equation ([5]) is the the main result of 
the present section, which will be used in the next section 
to analyze the existence of the lattice TDDFT. 

III. TDDFT ON A LATTICE 

The whole concept of TDDFT is based on the existence 
a one-to-one map between the time dependent density 



and the external potential. In this section we will prove 
the mapping and the w-representability theorems for the 
lattice TDDFT by adopting ideas recently proposed for 
the lattice version of TDCDFT^. 



A. Statement of the problem and the basic 
existence theorem 

In general we follow the NLSE approach to TDDFT- 
type theorie o 13 ' 17 ! 19 ' 20 . Let us first review the basics of 
this formalism in a more familiar continuum case and 
then discuss its modifications for many-body dynamics 
on a lattice. 

In a continuum system the starting point is the usual 
A^-particle Schrodinger equation 

id t \i,(t)) = [f + w + v{t)Mt)}, (ii) 

where T, W and v(t) are the kinetic energy operator, 
the interaction Hamiltonian and the external potential 
respectively. For a given initial state \ipo) the solution 
of Eq. (TTT1) determines the wave function \i/j(t)) and the 
density n(r; t) — (tjj(t)\h(r)\^(t)) as unique functional of 
the external potential, i. e. IV'MW) an d n{v]{t). Hence 
within the standard standard statement of the problem 
in quantum mechanics the linear Schrodinger equation 
defines a unique "direct map" from the external poten- 
tial and the initial state to the time-dependent wave 
function and the density: {v(t),\ipo)} h-> {n(t),\ip(t)}}. 
TDDFT assumes that in the above map the potential 
v(t) and the density n(t) can be interchanged. In other 
words, TDDFT relies on the existence of an "inverse 
map" from a time-dependent density and the initial state 
to the time-dependent wave function and a potential 
that produces the prescribed density: {n(t), \ipo)} ^ 
{v(t),\ip(t)}}. Constructively the inverse mapping can 
be defined as follows. Given the initial state |^o) and 
the density n(r; t) one finds a time-dependent wave func- 
tion by solving the Schrodinger equation ([lip , where the 
potential v(r;t) is not fixed externally, but determined 
self-consistently from the force balance equation 

iv.[n(r;t)V«(r;t)]=n(r;i) + g[V(*)](r). (12) 

In this equation the quantity q[tjj(t)](r) (the divergence 
of the stress force) is the following instantaneous bilinear 
functional of \i>(t)) 

qim}(r) = v(m\m,T+w]\m) (13) 

where j(r) is the usual operator of the current density. 
Apparently the initial state should be consistent with 
behavior of the density around initial time to to satisfy 
the following relations 

(Vo|ri(r)|Vo} =n(r;t ) (14) 
-V^o|j(r)|Vo)=n(r;t ) (15) 
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In this framework the proof of existence of TDDFT 
reduces to proving the uniqueness and existence of solu- 
tions to the nonlinear many-body problem of Eqs. ([TT]) . 
(j!2[) supplemented with an initial condition which satis- 
fies the consistency relations of Eqs. (TT4)) and (fT5j). 

Strategically the solution of the outlined nonlinear 
problem contains two major steps: (i) inverting the 
Sturm-Liouville operator in the left hand side of Eq. (fT"2")) 
to find the potential as functional of a given density and 
the instantaneous state, v[n(t), \ip(t))](r); and (ii) insert- 
ing this potential into Eq. (fTTj) and solving the result- 
ing NLSE. In a continuum proving the corresponding 
existence theorems for either step is a highly nontrivial 
and currently unsolved problem, although a significant 
progress has been made recentl y 21 ' 22 ' 35 . Below we re- 
formulate the problems (i)-(ii) for the lattice many-body 
theory and show that in this case a rigorous proof of ex- 
istence can be given. 

On a lattice the construction of the inverse map 
{n(t), \ipo}} ^ {v{t), consists of solving the dis- 

crete Schrodinger equation ([1]), where the potential v(r; t) 
is determined self-consistently from the lattice force bal- 
ance equation of Eq. ((5J) which we rewrite as follows 

J2krAmMr;t)=n(r;t)+ q m)](r)- (16) 

Here the functional q[ip(t)](r) is defined by Eq. © and 
we introduced the following notation 



T r , r ,p(r,rO-<5 r , r '^T r ^p(r,r") 



The initial condition for the nonlinear problem of Eqs. (UJ 
and (|16p should satisfy the consistency conditions, which 
follow from the definition of the density Eq. © and the 
continuity equation ((4]) 

N ^ |V>o(r,r 2 ,...,r w )| 2 =n(r;t ), (18) 

-2Im£ T r , r ,p Q (r, r') = n(r; to). (19) 

Equations (UJ, ([TB]). ([T5]l. and (O are the lattice 
analogs of Eqs. ([IT]). (fHj) . ([lij). and (JTSJ). A dramatic 
simplification of the lattice theory comes from the fact 
that both the Hilbert space H and the space V of lattice- 
valued potentials become finite dimensional, with the di- 
mensions M N and M, respectively. In particular the 
lattice TV-body Schrodinger equation ((T|) corresponds to 
a system of M N ODE, while the force balance equation 
(ITBl) turns out to be a system of M algebraic equations. 
In fact, Eq. (|16l) can be conveniently rewritten in a ma- 
trix form 



(20) 



where K is a real symmetric MxM matrix with elements 
fc r r / of Eq. (fTT)) , and V and S are M-dimensional vectors 



with components 

v r = v(r) and s r (n, ip) = h(r) + q[ip](T), (21) 

respectively. The A-matrix in Eq. (I2TJ1) is a lattice ana- 
log of the Sturm-Liouville operator m _1 VnV in Eq. (|12l) . 
Hence on a lattice the step (i) in solving our nonlinear 
many-body problem reduces to a simple matrix inver- 
sion, which can be performed provided the matrix K is 
nondegenerate. At this point it is worth noting that be- 
cause of the gauge invariance K matrix of Eq. (fTT)) al- 
ways has at least one zero eigenvalue that corresponds 
to a space-constant eigenvector. Therefore if V is the 
M-dimensional space of lattice potentials v(r), then the 
invertibility/nondegeneracy of K should always refer to 
the invertibility in an M — 1-dimensional subspace of V, 
which is orthogonal to a constant vector «c( r ) = C. 
In more physical terms this means that the force bal- 
ance equation (fT2"j) determines the self-consistent poten- 
tial v[n,rp](r) only up to an arbitrary constant. 

Now we are in a position to formulate and to prove the 
basic existence and uniqueness theorem on the lattice 
TDDFT. All statements of the Theorem 1 below refer to 
the lattice A-body problem defined in Sec. II. 

Theorem 1. (existence of the lattice TDDFT) — As- 
sume that a given time-dependent density n(r; t) is non- 
negative on each lattice site, sums up to the number of 
particles A, and has a continuous second time derivative 
n(r; t). Let be a subset of the ./V-particle Hilbert space 
H where the matrix K(ip) of Eq. (TT71) has only one zero 
eigenvalue corresponding to a space-constant vector. If 
the initial state tpo € f2, and at time to the consistency 
conditions of Eqs. (|18p and (fT9"f arc fulfilled, then 

(i) There is a time interval around £o where the nonlin- 
ear many-body problem of Eqs. (fTJ, (|16[) has a unique 
solution that defines the wave function ip(t) and the po- 
tential v(t) as unique Junctionals of the density n(t) and 
the initial state ^oi 

(ii) The solution of item (i) is not global in time if and 
only if at some maximal existence time t* > to the bound- 
ary of f2 is reached. 

Proof — By the condition of the theorem belongs 
to Q where K(ip) has only one trivial zero eigenvalue. 
Hence there is a neighborhood of tpo, such that for all ip's 
from this neighborhood the matrix K(tp) can be inverted 
(in the M — 1-dimensional subspace of V, orthogonal to a 
constant). In other words, we can solve the force balance 
equation (f2TJ|) as V — K~ 1 S and express (up to a con- 
stant) the on-site potential in terms of the instantaneous 
wave function and the density 

»[n,^(r) =^^- r 1 ,(VK(n,V). (22) 
Substituting this potential into Eq. ([T]) we obtain the 
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following NLSE, 

JV 

id t ip(r u ...,r N ;t) = K-] r s r ij>(ri, ...,r N ;t) 

3 = 1 r 

N 

XX 7 '' * x ;t) 

2—1 Xi 

+ ^2w r ^ rj ip(r 1 ,...,r N ;t). (23) 

j>i 

Equation (|23l) supplemented with the initial condition 
of Eq. ^ constitutes a universal nonlinear many-body 
problem, which determines the wave function in terms of 
the density. Formally it corresponds to a Cauchy prob- 
lem for a system of N-h = M N ODE of the following 
structural form 

ij> = F(iP,t), fl>(to) = fh, (24) 

where ip is a A^-dimensional vector living in the Hilbert 
space % and the right hand side is a nonlinear function 
of the ip's components. 

The nonlinearity of F(ip, t) in Eq. ([2^]) comes from the 
dependence of the potential v[n,ip](r) in Eq. (|2"2")l on the 
wave function, which, in turn, is determined by the func- 
tions K' 1 ^) and S(h,ip). Both k r ^(ip) of Eq. (IT71) . 
and s r (ri,'0) of Eq. (1211) are linear in the density ma- 
trices, and thus bilinear in if) forms. Therefore the po- 
tential v[n,ip](r) and, as a consequence, the whole right 
hand side F(ip, t) in Eq. (|24l) are rational functions of the 
components of the wave function. Moreover the denomi- 
nator of these rational functions never turns into zero for 
all ip € f2, which implies that in the function F(ip,t) 
satisfies a uniform Lipshitz condition. An explicit time 
dependence of F(ip,t) is determined by the time depen- 
dence of n{t) that is continuous by the condition of the 
theorem. Thus we conclude that for all ip € f2 the right 
hand side F(ip, t) in Eq. (fMf is Lipshitz in ip and contin- 
uous in time. 

After identifying the subset fi with the domain of Lip- 
shitz continuity we can directly employ the standard re- 
sults of the theory of nonlinear ODE. Namely, if the ini- 
tial state tpo £ the Picard-Lindelof theorem^ guaran- 
tees the existence of a finite interval to — 5 < t < to + 5, 
with 5 > 0, where the initial value problem (|24l) has 
a unique solution. This solution defines a unique map 
{n(t),ipo} ^ { u (*)> VK*)} locally in time, in accordance 
with the statement (i) of the theorem. 

The extension theorems for nonlinear ODE imply that 
a local solution can not be extended beyond some max- 
imal existence time t* > tg only in two cases: first, at 
t — > t* the solution becomes unbounded or, second, at 
t — > t* it reaches the boundary of fi. In our case the so- 
lution is guarantied to be normalized and thus bounded. 
Therefore we are left only with the second possibility, 
which proves the statement (ii) and completes the proof 
of the theorem. 



B. Discussion and comments on the existence 
theorem 

1. Definition of the v-representabihty subset Q 

According to the Theorem 1, any sufficiently smooth 
density n(r; t) is w-representable, at least locally, if the 
dynamics starts inside the subset SI of the Hilbert space. 
In general to ensure that a state ip belongs to f2 we need 
to check the invertibility of matrix K{ip), which, though 
possible in principle, may become difficult in practice, 
especially for lattices with a large number of sites. Is it 
possible to formulate simpler, but possibly more restric- 
tive criteria, which would guarantee the validity of the 
lattice TDDFT. 

One simple necessary condition immediately follows 
from the form of the matrix elements k r y in Eq. (|17l) . 
The matrix K{ip) is nondegenerate only if a lattice state 
ip is connected in a sense that any two sites can be con- 
nected by a line composed of links with nonzero values 
of T r ^Rep(r, r'). Indeed, for a disconnected state K(ip) 
takes a block-diagonal form and new zero eigenvalues, 
corresponding to piecewise constant in space eigenvec- 
tors, appear. The number of such zero eigenvalues equals 
to the number of disconnected regions on a lattice. We 
emphasize that a purely geometric connectivity of the lat- 
tice does not automatically guaranties w-representability 
- two sites connected geometrically by a nonzero hop- 
ping matrix element T r y can be disconnected in the 
above sense if for a state ip the quantity Rep(r, r') van- 
ishes. An explicit example of such a disconnected (one- 
particle) state on a connected 4-site tight-binding cluster 
has been recently proposed in Ref. [26| to demonstrate 
a possible non-w-representability in the lattice TDDFT. 
The authors considered an excited state with nodes on 
two opposite corners of a square formed by four sites 
(see Fig. 1 in Ref. [26h . The two nodes effectively sepa- 
rate the system into two disconnected parts. Therefore 
the matrix K acquires an extra zero eigenvalue and the 
Theorem 1 does not apply if the dynamics starts from 
such a state. In fact, one can show that this particular 
state lies precisely at the border of the w-representability 
subset fl. 

Obviously the connectivity of the lattice state is only a 
necessary, but not a sufficient condition for ip to be in fi. 
The reason is that for a connected state the quantities 
Tr,r'Rep(r, r') for different links may have different signs 
which can be responsible for extra zero eigenvalues of K. 
Hence the simplest sufficient condition is a connectivity of 
the state and positivity (or negativity) of T r r /Rep(r, r') 
for all lattice links. In other words, a state ip € 17 if its K 
matrix is primitive and does not have a block-diagonal 
form. This condition is easy to check in practice, but it 
appears to be quite restrictive. A less restrictive criterion 
that in many cases can still be checked easily, is a positive 
(negative) definiteness of K(ip). In Sec. IV we will show 
that this is exactly the case for a many-body ground state 
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on a connected lattice. Namely, if ipo is a ground state, 
then K(?Pq) is negative definite and thus ipo £ tt, which 
implies the existence of the lattice TDDFT for a system 
evolving from its ground state. 

2. Boundness ofh(r;t) in the lattice TDDFT 

A specific feature of quantum dynamics on a lattice, 
which narrows a class of w-representable densities, is a 
boundness of the time derivative of the densit y 20 ' 25 ' 28 ' 37 . 
Since the hopping rate along a given link is fixed to be 
T r r /, a link current of Eq. ^ can not exceed a cer- 
tain maximal value, i. e. |J(r,r')| < | J max (r, r')|, where 
I Jmax\ can be estimate d 20 ! 24 ' 28 using the Cauchy-Schwarz 
inequality 

|Jma X (r,r')| = 2|T r , r , /0 (r ) r')| < 2\T r ^\y/n{r)n{r'). 

(25) 

A physical density should satisfy the continuity equation 
© which imposes a bound on its time derivative, 

|n(r)| <]T|J max (r,r')|. (26) 

r' 

On the first sight the Theorem 1 does not say anything 
about the boundness of n(r; t). Therefore it is instructive 
to see how it can be deduced from the conditions of the 
theorem. First of all we note that if a solution of the 
universal NLSE exists, then the continuity equation is 
necessarily satisfied, which can be true only if our given 
density does not violate the bound of Eq. (|26|) . Accord- 
ing to the assumptions of the theorem the density should 
satisfy the consistency conditions, Eqs. (IT8l) and (|T9|) . 
and its second time derivative n(r:t) should be continu- 
ous in time for all t > t$. By imposing the condition of 
Eq. ([TTJf we explicitly require the boundness of n(r; to) at 
the initial time to, while the continuity of n(r; t) ensures 
that the physical bound of Eq. (|2"6"]l can not be violated 
immediately It is also worth noting that the bound- 
ness of n(r;t) is closely related to the invertibility of K 
matrix or, more precisely, to the connectivity of the in- 
stantaneous state tp(t). Indeed, the link current J(r, r') 
of Eq. ^ and the off-diagonal element fc rr < of Eq. (JT7J) 
are, respectively, the imaginary and the real parts of the 
quantity 2T r yp(r, r'; t). Therefore for any state i/j and 
r^r' the following identity holds true 

|J(r,r')| 2 + |fc r , r f = |J max (r,r')| 2 . (27) 

Equation (|2~?| shows that when the current J(r, r') 
reaches the maximal value of Eq. (1251) . fc r r ' turns into 
zero, which breaks the link between sites r and r'. Hence 
saturation of the bound in Eq. (f2l)f implies breaking all 
links attached to a site r. This site becomes disconnected 
from the rest of the lattice and the K matrix acquires an 
extra zero eigenvalue, indicating that the state tp is not 
anymore in the u-representability subset f2. Thus satu- 
ration of the bound on h(r) at some time t = t* auto- 
matically assumes that at this time the solution hits the 



boundary of f2. This behavior is in a clear agreement 
with the statement (ii) of Theorem 1. 

IV. EXPLICIT ILLUSTRATION: ONE 
PARTICLE ON A TWO-SITE LATTICE 

This section is aimed at illustrating the general NLSE 
scheme using a simple exactly solvable example - one 
particle on two sites. In spite of its simplicity this ex- 
ample contains practically all features of the most gen- 
eral JV-body case, and thus displays all subtle points of 
the general formulation in a clear and transparent form. 
Therefore it is advisable to read this section to get bet- 
ter feeling of the formalism presented in Sec. III. Another 
explicit example of NLSE for an interacting two-particle 
system is given in Appendix. 

Consider a particle living on a two-site lattice. The 
state of the system is described by the one-particle wave 
function ip r (t), where the coordinate r takes values 1 or 
2 corresponding to the two lattice sites. 

The dynamics of the system is described by Eq. (P) 
where the number of the particles N — 1 and there is no 
interaction term in the right hand side. Therefore Eq. (fTJ) 
reduces to the following system of two ODEs. 

idtiPx = uiVi - (28a) 
id t ip 2 = v 2 ip 2 -Tilj 1 , (28b) 

where T is a real hopping rate, and v\ and V2 are the 
time-dependent external potentials on sites 1 and 2, re- 
spectively. The system of Eqs. (f2"5)) determines the com- 
ponents of the wave function ip± and as functionals of 
the external potential and the initial state ^1,2(0) (here 
we set to = 0). 

To find the wave function as a functional of the density 
we have to construct the proper NLSE, and for this one 
needs an additional equation which relates the potential 
to the density and the wave function. In the general 
framework of Sec. Ill the force balance equation of Eq. ([5]) 
[or, equivalently, Eq. (|20|) ] serves exactly for this purpose. 
For two sites K is a 2 x 2 matrix and therefore Eq. (j2"0|) 
takes the form 

f? 12 k r)( vi ) = ( hl t qi ) (29) 

where k rr > and q r can easily be derived from the general 
definitions of Eqs. (fT?)) and @, respectively, 

fc 12 = 2TRep 12 , (30) 
qi = -q 2 = -2T 2 {\^ 2 \ 2 -\^\ 2 ). (31) 

The link density matrix p\2, which is in general deter- 
mined by Eq. ([5]), in the present one-particle case re- 
duces to a simple product p\2 = ipi^2- The 2x2 matrix 
in right hand side of Eq. ([29]) is the K matrix extensively 
discussed in the previous section. Obviously it always has 
a zero eigenvalue corresponding to a space-constant po- 
tential v\ — V2 = C. If there is no other zero eigenvalue 
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we can invert K in a space perpendicular to the constant 
vector «c(r) = C. Being perpendicular to Uc(r) simply 
means that on-site potentials sum up to zero, which for 
two sites implies 



-v 2 = v. 



(32) 



This equation defines a f-dimensional subspace of V 
where K can in principle be inverted. In the present two- 
site case the K matrix is invertible if fci 2 7^ 0. Therefore 
the v-represenatability subset of the Hilbert space is 
defined by the following simple condition 



Re[^ 2 ] ^ 0. 



(33) 



For all states satisfying the condition Eq. (f3"3"|) we can 
invert K matrix in Eq. ([29|i and find the potential v — 
v\ = —V2 as a functional of the density and the wave 
function, 



n 1 -2T 2 (|^ 2 | 2 -|^ 1 | 2 ) 
4TRepi 2 



(34) 



where we substituted explicit expressions for fci 2 and qi 
from Eqs. (00) and (|3T|). 

The final NLSE is obtained by inserting the potential 
of Eq. (|34) into the Schrodinger equation ([28)) 



n 1 -2T 2 (|^ 2 | 2 -|V; 1 | 2 ) 
Z ^ = iTRe^H * ~ ^ (35a) 

4TRe[V> 1 *?/' 2 ] ^ ^ v 7 

This system of equations perfectly illustrates the generic 
structure of NLSE appearing in the TDDFT context. 
Firstly, as described in Sec. Ill, the nonlinearity is always 
a rational function with enumerator and denominator be- 
ing bilinear forms in the components of the wave function 
-0. For all ip € f2 [i. e. for ip satisfying Eq. (|3"3")) ] the de- 
nominator never turns into zero. Secondly, the explicit 
time dependence enters NLSE only via the second time 
derivative of the density n r (t) that is assumed to be con- 
tinuous. The above two properties ensure that the right 
hand side of our NLSE is Lipshitz in tp and continuous 
in t. By the Picard-Lindelof theorem this guarantees the 
existence of a unique solution to Eqs. (1351) for any initial 
state ?/>(0) from il. 

However, this is not yet the whole story. Since the den- 
sity enters the equations only via h r (t), our unique so- 
lution to NLSE, in general, will reproduce correctly only 
the second time derivative of the prescribed density. The 
whole externally given density n r (t) is recovered from 
NLSE if the dynamics starts from a special manifold of 
the "density-consistent" initial states which are defined 
by the consistency conditions of Eqs. (fT5)l and (flU)) . 

To proceed further with our example we represent the 
wave function in the polar form 

Mt) = |^iWI^ (t)/2 , Mt) = \Mt)\e~ Mt)/2 , (36) 



and substitute it into the consistency conditions. As a 
result Eqs. (fT5]) and ([H?|). respectively, simplify as follows 



|<M0)| = r = {l,2}, 

hi(0) = 2iy«i(0)n 2 (0) sin 95(0), 



(37) 
(38) 



where n r (0) is the (prescribed) initial density, and n r (0) 
its initial time derivative. The first condition, Eq. (|37l) . 
uniquely determines the modulus of the allowed ini- 
tial states in terms of the initial density. Finding the 
"density-consistent " initial phases is a bit more tricky as 
the right hand of Eq. (I3"51) is not a single valued function. 
Equation (|38l) has two solutions which can be written in 
the following form 



tp^ (0) = arcsin | 
^-)(0)=7r-^ , 



ni(0) 



27yn 1 (0)n 2 (0) 



(39) 
(40) 



where arcsin stands for the principal value of the inverse 
of sin. In other words, </>o defined after Eq. (|3"9")l is a so- 
lution to Eq. (|38|) in the interval [— f , f ]■ The existence 
of two solution to the consistency conditions means that 
for a given density in our simple system the manifold of 
density-consistent initial states consists of the two follow- 
ing wave functions 

^ +) (0) = vWMe^, 4 +) (0) = v/^(0) e -%l) 

^- } (0) = v^Oje-*^, 4~ } (0) = -v^(0y%2) 

where we disregarded an irrelevant common phase factor. 
By substituting Eqs. (l4"TT) and (|4"2")l into the condition 
of Eq. (|33]) we find that the initial states t/i (±) (0) S fl 
provided 4>o ^ ±ir/2. Obviously this puts a restriction 
on the initial values of n and h. Equation (131?)) tells us 
that the condition </>o ^ ±7r/2 actually ensures that n(0) 
is properly bounded [see Eqs. (|25|) and ([26])]. 

It is interesting to note that if ri(0) = 0, than <fio = 
and the density-consistent initial states Eqs. (|41l) and 
(1421) can be viewed as the ground (symmetric) and the 
excited (antisymmetric) states of a dimer in the presence 
of some static potential. 

Now we can solve NLSE of Eq. (|35p starting from one 
of the allowed initial states. The solution should return 
the wave function and the potential as unique functionals 
of the given density n r (t). Inserting the polar representa- 
tion Eq. (|3l))) into Eq. (|3"5"j) we observe that the following 
form of ipr(t), 

Mt) = AM^w/ 2 , Mt) = VMtY~ Mt),2 i (43) 

solves NLSE if the time-dependent phase ip (t) satisfies 
the equation 



hi(t) = 2T '\J 'ni(i)n 2 (£) smtp{t). 



(44) 



For each initial state from the set of Eqs. (|4T|) and (|42l 
this equation has a unique solution provided, the condi 
tion fci 2 (V0 ^ is fulfilled. 
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Assume that we started from the state ip\ 
Eq. (0TJ). Then the solution to Eq. @J reads 



(+) 



(0), 



<P 



\2Ty/nTir 2 



(45) 



This equation together with Eq. (|43|) gives the wave func- 
tion as a functional of the density. As long as this solution 
exists, the element k\2 = 2TRepi2 of the K matrix stays 
positive 



■„(+) 

»12 



1T^Jri\ni cos </? = W4T 2 niri2 



(46) 



To find the potential v(i) as a functional of the density 
we insert Eqs. fl33J and dMl into Eq. (|33]). The results 
takes the following form 



fii - 2T 2 (m - na) 



111712 — TK 



(47) 



This functional reproduces with the result obtained in 
Ref. In addition there is another solution that corre- 
sponds to another density-consistent initial state. 

If we start from the second initial state, ip[ \o) of 
Eq. (|42|) , we should take the second solution of Eq. (|44)) 
for the phase, namely 



. f ni 

(p = it — arcsm 



In this case k\2 changes a sign, 



u 12 



4T 2 niri2 



(48) 



(49) 



which implies that the sign of the potential v is also re- 
versed 



(-)r l ni 



2T 2 ( ni -n 2 ) 



n\U2 — n\ 



(50) 



Thus, for different initial conditions the NLSE machinery 
produces unambiguously different functionals v[n] and 
if>[n]. This nicely displays the initial state dependence 
in TDDF T 38 ' 39 , although in the present case the depen- 
dence is very simple. 

For each density-consistent initial state the unique so- 
lution to NLSE of Eq. (1551) exists as long as the bound- 
ary of the subset fi is not reached. This happens if 
fci2(V0 = 0, i. e. when the expression under the square 
root in Eqs. (l46l) or (l49l) turns into zero. In agreement 
with the general discussion in Sec. Ill, at this point the 
bound on the time derivative of the density, 



| Al| < 2Ty/nin 2 , 



(51) 



is saturated. 

In our simple model we can also visualize and 
completely characterize the geometry of the v- 
representability subset fi in the Hilbert space H. Since 
we have an effective 2-level system the projective Hilbert 




FIG. 1. Each normalized state in the Hilbert space % maps 
to a point on the Bloch sphere. The north |1) and the south 
|2) poles correspond to the particle on sites 1 and 2. The line 
k\2 = divides the sphere into two (left and right) hemi- 
spheres corresponding to two disconnected parts of the v- 
representability subset £1. 



space can be represented by a 2-sphere, known as a Bloch 
sphere. Specifically, after taking out a common phase fac- 
tor, each normalized state from % is mapped to a point 
on a 2-sphere in R 3 (see Fig. 1) by parametrizing the 
wave function as follows 



cos6»/2< 



iip/2 



1) + sin 9/2 



,-i(p/2 



|2), 



(52) 



where |1) and |2) are the orthogonal states correspond- 
ing to the particle residing on sites 1 and 2, respectively. 
In this mapping the moduli of the on-site amplitudes are 
represented by the azimuthal angle 9, while the phase dif- 
ference ip corresponds to the polar angle in spherical co- 
ordinates. As we can see from Fig. 1 the line fci2(V0 — 
divides the projective Hilbert space into two, left and 
right, hemispheres with fci2 > and fci2 < 0. The two 
hemispheres represent two disconnected parts of the v- 
representability subset ft, separated by the boundary line 
ki2 = 0. The boundary line contains all states for which 
the statements of the Theorem 1 do not hold. Starting 
from any point on those hemispheres, i. e. from a state 
ip £ f2, we uniquely recover the time evolution of the sys- 
tem with a given density by solving NLSE of Eq. (|35l) . 
As long as the trajectory stays within the original hemi- 
sphere and does not touch the boundary, the one-to-one 
density-to-potential map exists with the functional v[n] 
given by Eq. (|4T|) or by Eq. (|50")l , depending on the hemi- 
sphere. Whether it is possible to construct a unique and 
universal analytic continuation for crossing the bound- 
ary and covering the whole subset f2 is an interesting 
question, which can not be answered at the level of The- 
orem 1. 
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V. TIME-DEPENDENT u-REPRESENTABILITY 
FOR A SYSTEM EVOLVING FROM THE 
GROUND STATE 

In this section we return to the most general case, and 
show that the ground state of a lattice iV-particle system 
always belongs to the i;-representability subset f2. This 
implies that the lattice TDDFT is guaranteed to exist if 
the dynamics starts from the ground state. 

Assume that ip^ = \k) form a complete set of eigen- 
states for a lattice many-body Hamiltonian describing 
./V-particle system in the presence of a static scalar po- 
tential i>o(r). Let Vo = |0) be the ground state. We are 
going to prove that the matrix K(ip ) evaluated at the 
ground state is strictly negative definite in the subspace 
of potentials that are orthogonal to a space-constant vec- 
tor V c , i- e., 

V T K(ip )V = v(r)k ry v(r') < 0, (53) 

r.r' 

for all M-dimensional vectors V — {v(r)} which satisfy 
the following orthogonality relation 

yTy c =C$>( r ) =0, (54) 

r 

where V T stands for a transposed vector. 

Our starting point is the /-sum rule (see, e. g., Ref.l40l) 
for the density-density response function Xr,r' (w) = 
((n r ;nr')) u : 

2 poo _ 

/ wlmXr r ' (oo)duj = i(0\ [n r , n r >] |0). (55) 

7T Jo 

To calculate the commutator in the right hand side of 
Eq. (|55p we switch to the second quantized representation 
and write the equation of motion for the density operator 
n r = ajd r , 

= i ^2(T ry ala r < - T r ', r a^a r ), (56) 

r' 

where a r and aj are the on-site annihilation and creation 
operators. Equation (|56[) is nothing but the operator 
form of the continuity equation. Using Eq. (|56[) one can 
easily calculate the commutator entering the right hand 
side of Eq. ([55]). 

i[h r ,n r >] = — T r yala T > + S rr > T r , r "a£a r " + h.c. (57) 

r" 

Taking the ground state expectation value of this equa- 
tion and comparing the result with Eq. (|17[) we find 

i(0\[fi T ,h r ,}\0) = -k ry . (58) 

Therefore the right hand side of the lattice /-sum rule 
is identified with the K matrix entering the definition of 
the u-representability subset 4 ^. 



On the other hand, for the imaginary part of den- 
sity response function we have the following spectral 
representation 40 

oo 

J u)lmxr,r'(u)duj = — 27rRc a;fco(0|n r |fc) (fc|n r /|0), 
o k 

where uiko = Ek — Eq are excitation energies of the sys- 
tem. Substitution of Eqs. J5SJ) and ^ into Eq. flSSJ) 
leads to the spectral representation for the elements of 
K matrix 

kr, T ' = -4:ReJ2uj ko (0\n r \k)(k\n r ,\0). (60) 

k 

Finally, inserting k r r > of Eq. (f60|) into the left hand side 
of Eq. (I53|) . we arrive at the following remarkable result 

V T k{^ )V = -4 Wfeo | J2 «(r)(0|n r |fc) J' 

k r 

= -4^cj fco K0|fi|fc)| 2 <0 ! (61) 

k 

where v is a many-body operator corresponding to the 
potential v(r), 

v — v(r)n T . (62) 

r 

Let us show that the equality in Eq. (|6T|) holds only for 
a space-constant potential uc( r ) = C Since each term 
in in the sum in Eq. (|6Tj) is non-negative, the result of 
summation is zero if and only if 

(0|u|fc) = 0, for all k ^ 0. (63) 

Physically the right hand side of Eq. (|6"Tj) is proportional 
to the energy absorbed by a system after a small ampli- 
tude pulse of the form v(r; t) — v(r)5(t). Then the condi- 
tion Eq. (|63|) simply states that nothing is absorbed only 
if the potential v(r) does not couple the ground state to 
any of the excited states. 

Assume that Eq. (|63[) is fulfilled and expand a vector 
v\0) in a complete set of states {\k)} 

£|0)=]T|fc)(fcH0)H0)(0|{;|0)^A|0). (64) 

k 

Therefore the condition of Eq. (|63| implies that the 
ground state |0) is an eigenfunction of the operator v. 
Since v corresponds to a local multiplicative one-particle 
potential this can happen only if the potential is a con- 
stant. For clarity we write Eq. (|64| in the coordinate 
representation 

N 

^2v(Tj)ipo(n, ...,rjv) = AV>o(ri,...,rjv). (65) 

3=1 

Obviously this equation can be fulfilled only if the func- 
tion v(r) takes the same value X/N on all lattice sites, 
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which corresponds to a space-constant potential. A no- 
table exception is a geometrically disconnected lattice 
consisting of a several pieces that can not be connected 
by a path composed of links with nonzero T r r i . In this 
case the arguments of the wave function ^>o( r i> r Jv) 
form "disconnected groups" of coordinates corresponding 
to particles residing in disconnected parts of the system. 
The coordinates of different disconnected groups take val- 
ues in "non-overlapping" parts of the lattice. Since the 
number of particles in each part (number of coordinates 
in each group) is fixed, Eq. ([551 can also be satisfied with 
a piecewise constant potential. 

Therefore we arrive at the following conclusion: for a 
connected lattice Eq. (f63|) is fulfilled, and the inequality 
in Eq. (|6ip is saturated only for a constant in space po- 
tential. For all potentials which are orthogonal to a con- 
stant in a sense of Eq. ([54")) the strict inequality in Eq. (f6Tj) 
takes place. This means that matrix K(ipo) is negative 
definite and thus invertible in the M — 1-dimensional sub- 
space of V orthogonal to a constant vector Vc- In other 
words, the ground state of TV-particle system on a con- 
nected lattice does belong to the u-representability subset 
fi. This result combined with the general existence theo- 
rem of Sec. Ill proves the following particular version of 
the time-dependent w-representability theorem. 

Theorem 2. - Let the initial state i/'o for a time- 
dependent many-body problem on a connected lattice 
corresponds to a ground state in the presence of some 
scalar potential Vo(r). Then any density n(r;t) which 
satisfies the consistency conditions of Eqs. (fTS]) and (TT9]) 
and has a continuous second time derivative is locally 
f-representable. 

It is worth noting that the above ground state based 
argumentation can be straightforwardly extended to a 
thermal equilibrium state. In this case Eq. (|61j) takes 
the form 

k>l 

where (3 is inverse temperature, K@ is the K matrix eval- 
uated for the thermal equilibrium state, Z is partition 
function, and LJki = Ek — E%. By the same token the 
form defined by Eq. (|66p is strictly negative for all po- 
tentials orthogonal to a constant vector. Therefore the 
Theorem 2 should also apply to the ensemble version of 
TDDFT based on the von Neuman equation for the N- 
body density matrix. Of course in this case we need to 
prove an ensemble extension of Theorem 1 , but currently 
this also looks relatively straightforward. 

VI. CONCLUSION 

In conclusion we presented a rigorous formulation of 
the lattice TDDFT. On the technical side we mapped 
the problem of existence of TDDFT to the Cauchy prob- 
lem for a special NLSE and proved the corresponding 



existence and uniqueness theorem using standard results 
from the theory of nonlinear ODE. As one could expect 
from the previous works2£-^2ML the lattice TDDFT is 
not unconditionally valid. Similarly to the lattice TD- 
CDFT, apart from a certain quite mild restrictions on 
the time dependence of the density, there are also limi- 
tations on possible many-body initial states. Namely, al- 
lowed initial states form a special w-representability sub- 
set ft of the Hilbert space H. Most physical states should 
be in ft. In particular we demonstrated that this is true 
for any ground state on a connected lattice. Therefore 
in practice for most non-exotic physical situations the 
lattice TDDFT (and TDCDFT) should be valid. In this 
context it is interesting to note a recently observed exam- 
ple of a violation of the non-interacting v-representability 
for a strongly biased interacting resonant level model^. 

Our proof of the basic theorem for the lattice TDDFT 
as well as a recent proof for the lattice TDCDFT— essen- 
tially rely on the concept of the u-representability sub- 
set £1 in the Hilbert space. It would be very interest- 
ing and useful for the future to carefully study and to 
characterize the structure/geometry of CI. Formally the 
v-rcprcscntability subset is defined as a part of % where 
detf-RT] 5^ 0. Similarly to the explicit example of Sec. IV, 
in general ft should consist of two regions with differ- 
ent signs of det[AT], which are separated by a surface 
det[AT] = 0. What is the geometry of each part for a 
general quantum system? Are they simply- or multiply- 
connected? Is it possible for system driven a physical 
potential to cross the surface detfif] = or it is forever 
confined to one initially fixed subregion of f2. Answering 
these questions will definitely deepen our understanding 
of TDDFT-related theories as well as quantum dynamics 
in general. 
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Appendix: Universal NLSE for two interacting 
particles on a two-site lattice (Hubbard dimer) 

Below we present the derivation of the universal NLSE 
for a system of two interacting spin- 1/2 fermions on a 
two-site lattice. This gives another explicit example il- 
lustrating a general structure of nonlinear equations ap- 
pearing in Sec. III. 

Since the dynamics of a triplet state on two-sites is 
trivial we concentrate on a singlet state. The spatial part 
of the singlet 2-particle wave function is symmetric with 
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respect to the permutation of coordinates. It has three 
components, ipu, ^22, and ^12 = 1P21, describing differ- 
ent distribution of the particles over the sites. Because 
only the potentials orthogonal to a constant are relevant 
(see Sec. IV) we assume from the very beginning that 

V\ = —V2 = V 

The Schrodinger equation Eq. ((T|) for two particles in 
a singlet state on a two-site lattice takes the form 

idtipn = —Tipiz + (2v + ^11)^11 (A. la) 

idtipi2 = —Tipu + W121P12 - Tip22, (A.lb) 

id t ^p22 — —Ti>\2 + (—2V + W22)V"22- (A.lc) 

To construct the NLSE for this system we need to 
substitute the potential from the force balance equation 
which relates on-site potential to the density and the 
wave function. As we have a two-site lattice the struc- 
ture of the force balance equation coincides with that of 
Eq. derived in Sec. IV. We only have to substitute 
into Eq. ([29]) the values of k i2 = 2TRepi 2 and q r calcu- 
lated for the two particle system [by means of Eqs. ([5]) 
and ©] 

k V2 = 2TRe Pl2 = 4TRe[^ 11 V12 + ^22] , (A.2) 



Ql =- 4T 2 (|V 22 | 2 - |<M 2 ) + 4T(Re[V>^i2]Ki - w 12 ) 
- Re[-02 2 V'i2]('!«22 - W12) ) = -92- (A.3) 



the potential v in terms of density and the wave function 



As we know K matrix has a trivial zero eigenvalue 
corresponding to the constant potential which is already 
projected out by setting v — v\ — —v 2 ■ If 2TRepi 2 7^ 
we can invert K matrix in Eq. (|29[) . which determines 



v(n,ip) 



~4TRepi 2 ' 



(A.4) 



Finally by substituting v of Eq. (|A.4|) into Eq. (|A.1[) 
we arrive at the universal NLSE for this system 



ni + qi 



+ wiijipn 



2TRe[^ 1 Vi2 + ^1*2^22] 

- TVi2, (A.5a) 

idtipi2 = w\2^i2 - Tipn - T1P22, (A. 5b) 
n\ + qi 



idti>22 = ( 



2TRe[V> 11 V'i2 + ^22. 



+ W 2 2 S j'>p22 



(A.5c) 



We again explicitly see a system of nonlinear ODEs with 
the nonlinearity of the rational form. As one expects 
on the general grounds (see Sec III) inclusion of inter- 
actions does not introduce any conceptual modification 
in comparison with the simplest one-particle case con- 
sidered in Sec. IV. In the subset ft of the Hilbert space 
where Repi 2 ^ all the terms in the right hand side of 
Eqs. (|A.5al) and (|A.5c[) are infinitely differentiable with 
respect to components ip rr > of the wave function. Hence 
the whole right hand side stays Lipshitz continuous, pro- 
vided if] G 0. If the density has a continuous second time 
derivative n, and the initial state ip(0) € ft is density- 
consistent in a sense of Eqs. (fll^ - pUl) . then there is a 
unique solution to NLSE of Eqs. (|A.5[) . Therefore there 
exists a one-to-one mapping between the density and the 
potential in the original two-body problem. 
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